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Correlator Compensation Requirements for Passive 
Time-Delay Estimation with Moving Source 

or Receivers 

WILLIAM B. ADAMS, JOHN P. KUHN, and WILLIAM P. WHYLAND, member, ieee 



Abstract- An analysis is given of the effects of source or receiver 
motions on the output of a cross correlator used to estimate the source 
time-delay difference across a baseline. For both narrow-band and 
wide-band correlators, the need for time-varying correlator delay com- 
pensation is quantified for the case of a quadratic signal delay-difference 
variation. Two useful concepts which emerge are 1) the three-dimen- 
sional delay/delay-rate/delay-acceleration mean ambiguity function of 
the source signal, and 2) the equivalent r- domain filter, which trans- 
forms the source autocorrelation function into the mean output of the 
mismatched correlator. 

The required correlator compensation is approximately a quadratic 
delay modulation matched to the input delay-difference function. For 
3 dB peak output loss with a narrow-band signal, die maximum allow- 
able delay-rate mismatch will produce 158° of phase rotation at the 
centroid frequency f 0 during the correlation integration time 7\ while 
the maximum allowable delay-acceleration mismatch will produce 156° 
of quadratic phase rotation at fo between the center (f = 0) and each 
edge (t = ±T/2) of the correlator integration window. For broad-band 
signals, the mismatch tolerances become about 11 percent tighter. 



I. Introduction 

GIVEN two sensors separated by a baseline, passive esti- 
mation of the time-delay difference between signals re- 
ceived from a common source may be accomplished by cross 
correlating the appropriately filtered [1] sensor outputs. 
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When source or receiver motions cause the delay difference to 
vary, during the correlator integration time T, by more than 
the correlation time width of the source signal, then the out- 
put correlogram is degraded (smeared) unless the correlator 
implements a compensating delay modulation during T. In 
many cases the simple expedient of reducing T causes an un- 
acceptable reduction in processing gain. 

For a linear delay-difference variation during T, the opti- 
mum correlator delay modulation is also linear in time during 
T, so that one receiver output is time scaled (compressed or 
expanded) by a constant factor before cross correlating [2] . 
For a narrow-band signal, this compensation can be approxi- 
mated by a constant frequency shift [3], [6], [7] . 

For a quadratic delay-difference variation during T, the 
optimum correlator delay modulation is approximately 
quadratic in time during T, so that one receiver output is time 
scaled by a linearly varying factor before cross correlating. 
For a narrow-band signal, this compensation can be approxi- 
mated by a linearly varying frequency shift [3] . If the corre- 
lator can only implement a linear delay modulation (or its 
narrow-band approximation by a constant frequency shift), 
then the maximum allowable value of Tfor a specified correlo- 
gram degradation becomes a function of the geometry and 
motion scenario, together with the centroid and width of the 
source spectrum [3], [4]. 

All these effects for either broad-band or narrow-band corre- 
lators are treated in a unified manner by employing analytic- 
signal representation, assuming a stochastic source emission, 
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and examining the mean (ensemble average) correlogram out- 
put as a function of errors (mismatches) in the correlator 
delay-modulation schedule. Correlator compensation require- 
ments, for a specified degradation of 3 dB in the mean peak 
correlator output, are then obtained as error budgets on the 
allowable mismatches in compensating for delay-diffeience 
rates and accelerations (which at any reference frequency, and 
in particular at the centroid frequency f 0 of the source spec- 
trum, may be translated into allowable values of frequency 
difference and frequency-difference rate). 

n. Mean Correlator Output with Quadratic 
Delay Variations 

A, General Approach 

Assume the source emits a waveform s(r), which is a sample 
function of a zero-mean stationary random process with auto- 
correlation function (acf ) 

«,(r)=^{s(f)sa + r)}. (1) 

Neglecting medium distortion and attenuation, the signal 
received at sensor (*= 1,2) is delay-modulated because 
of source or receiver motion 



fft(3 -«['-*(;)]. 



P) 



If the propagation delays r k (t) are assumed to vary quadrati- 
cally with time during the observation interval 7 0 bs> 

r k (t) = r ok + r' ok t + r" ok t */2, 1 f | < r abs /2, (3) 

then the delay difference is also quadratic in time 

A(0 = Tj(0-r 2 (0 

= Ao+A^ + ASr 2 /2 (4) 

where 

A 0 = T 0 i ~ T <n 
A^ = Tor - t'm 

Ao*=r 0 W«. ( 5 > 
The cross correlator must then implement a time-varying 
delay-modulation schedule during its integration time T to 
compensate for the moving correlation point. We assume a 
quadratic 1 correlator delay-modulation schedule given by 

Ac(r) = r c + V + *c*V2 9 1 1\ < T/2 (6) 

so that when the correlator inputs, including additive noises, 
are 

y{t) = s[t-r 2 (f)] + « 2 (f), 
then the correlator output is 



The correlator attempts to compensate for Ai and A£ by a 
proper choice of \. and a Cy after which it performs (e.g., by 
parallel processing) a scan in r c to find the delay at which 
F is a maximum. The value is its estimate of the delay 
difference A 0 between the received signals at the instant t - 0 
(taken as the center of its integration window). 

Since F is a random function of r c because s> n x and n 2 are 
sample functions of random processes (assumed independent), 
we require a large product (T • A^nai) so that the time average 
performed by the correlator gives a useful approximation to 
the ensemble average over the noise and signal distributions, 
i.e., so that F versus r c is a good sampling image of the signal 
acf with a delay shift. The correlator compensation require- 
ments may then be determined by studying the mean (en- 
semble average) correlator output [5], [6] R a = E Stni9Hi {F} 
as a function of the mismatches in choosing t c , X c , and oc c . 
Combining (7) and (8) and averaging over n x and n 2 yields 
R d =E s {F\n l =n 2 = 0}, which expresses the well-known [5] , 
[1 1 ] result that additive zero-mean independent noises do not 
affect the mean cross-correlator output. Thus for our pur- 
poses the additive noise will be ignored. 

B. Analytic-Signal (Preenvelope) Notation 

In discussing both wide-band and narrow-band correlators 
it is convenient to use analytic-signal or preenvelope notation 
[6], [8], [9], where O denotes [()+/(")] with C) the 
Hilbert transform of ( ). 

Thus the real source waveform is 



(7) 



i r m 

F=F{r c , K, «c) = V *(0 yl* - ^(01 dt. 

1 J -T(2 



(8) 



s(f) = Re[?(f)] (9) 
where 

s(f)=»(0 +/*(') 0°) 
and the complex acf of s{t) may be shown to be [6] , [8] , [9] 

Kj<r) = *?{s*(')s(' + r)} 

= 2/?,(T)=2[*,(r) + /£ 4 (r)] (11) 

which Is twice the preenvelope of R t (j). 
With a source mean power density spectrum W„(f), then 



*,(T)' 

and [8] 



11 * 



<J)*™df 



= 2 f W 3 



(f)e /2 " fT df. 



(12) 



(13) 



1 This delay schedule can give perfect compensation for a constant or 
linearly-varying delay difference during T, but can only give a good ap- 
proximation to perfect compensation when rJk * 0 (Appendix B). 



C Mean Output of Delay-Modulated Correlator 

Using analytic-signal notation, the basic expression for the 
mean correlator output R d is now derived. Referring to the 
functional diagram shown in Fig. 1, the analytic source signal 
s(f) experiences the time-varying propagation delays T k (t) 9 
k = 1,2, and appears at correlator input port as s k (f) = 
s(t - Tjt(r)). The correlator imposes the compensating delay 
modulation A c (f) on s 2 (t) to yield s t [t- A c (f)] = 

s{t- A c (r)- r 2 [t- A c (r)]}» and then evaluates the time 
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Fig. 1. Schematic for mean output of delay -modulated correlator. 



average over (-772, T/2\ of the product Si(f)s?(0- The 
mean correlator output R d is then half the ensemble average 
(over s) of this finite time average, which from (1 1 ) is seen to 
be 



r 2r J, 



T72 



-Tf2 

-772 

2r 



1 f 772 



-H 



7/2 
r/2 



^[e(0] * 



(14) 



'-r/2 



where 

e(/) = A c (t) - {7,(0- r 2 [f - A c (f)] } . (15) 

is the negative of the delay difference at the input to the 
multiplier. 

Evaluating (15) with (3) and (6), we obtain after some minor 
approximations (Appendix B) 

e(0 ~ t c - A 0 ) + (0 2 K - A(, ) t + (&<*, - A n 0 ) f 2 /2 



= fe(T + X^ + af 2 /2) 
where we define 

02 = 1 - T02 
T = T C - A 0 /fe 

X=X C - A^/fla 
a = a c - AS/fe 
so that (14) is 



(16) 



(H) 



~ - 1 r 772 - 

K*=*<*0-, A,a) = — ^[Wr + Xf+o/ 2 /2)]^ 
7 J -r/2 



(18) 

which may be interpreted [7] as a three-dimensional mean 
ambiguity function of the source waveform. For our purposes 
it is simply the mean preenvelope of the correlator output, as 
a function of correlator mismatches in delay (t), delay rate 
(X), and delay acceleration (a). With perfect compensation in 
delay rate and delay acceleration, X = a = 0 so that 

£ d (T,0,0) = £,(feT). (19) 



For a wide-band correlator, the mean real output is given by 
the real part of Rd(?> X, a) from (18). For a narrow-band cor- 
relator this output versus r has the form of a modulated carrier, 
and since the delay as measured by the "carrier" phase is too 
ambiguous to be useful, we rely instead on the modulated 
source bandwidth for delay determination [7] by finding that 
T cm which maximizes the envelope \Rd(r, X, a)\. 

D. Equivalent Filter for Mismatched Correlator 

A useful concept in understanding the effects of correlator 
mismatches on the mean output R d is that of the equivalent 
filter [5]. If we put (13) in (18), interchange the order of 
integration, and change variable 

ft/-*/. (20) 
then we obtain 



where 



1 f r/2 

1 J -TJ2 



(21) 



(22) 



The bracketed term in (21) is the Fourier transform of 
^y(02 r )> and H d (f X, at) is an equivalent filter in the follow- 
ing sense: with r analogous to a time variable, then R d (r, X, a) 
considered as a function of r is the result of passing the "time" 
waveform R 5 (P 2 r) through a filter with frequency response 
Hdif> \ a ) r-domain impulse response 

Mr, K <*) = f H d (/, X, a) e^ df (23) 
J -00 

Note that H d is not a physical filter, but rather a mathe- 
matical analogy, as indicated by Fig. 2. 
By completing the square in the exponent in (22), we obtain 

e ~MxV/a) f 17 2X\ 1 



2jc 



+ Z 



(24) 



where 

x = Ty/fal2 
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Fig. 2. Equivalent filter for delay-modulated correlator. 



and 



JrU 
o 

is the complex Fresnel integral, obeying 
Z(-a)=-Z( W ) 
Z(j«)«/Z*(a). 
For special cases we obtain directly from (22), 
" d (/,0,0)=l 
tf d (/,A,0) = sinc (X/r), 
sin7rw 



(25) 



sine (w) = 



mi 



while from (24) we obtain 
#<,(/, 0,a)=^. 



(26) 
(27) 

(28) 
(29) 



III. Correlator Compensation Requirements for 
Delay-Difference Rates and Accelerations 

A. General Approach 

For both narrow-band and broad-band source spectra, we 
consider peak loss, resolution, and possible bias in the mean 
correlator output as a function of mismatches in the corre- 
lator delay-compensation parameters. Most of the results are 
well known when expressed in terms of range variation, rather 
than (as here) delay-variation coefficients. For physical inter- 
pretation and comparisons, we assume that the range from the 
source to sensor "k" (k= 1, 2) varies quadratically during the 
observation interval 

r*(0 = *>* + 4t* + 'S*'V2^ IH < W2. (35) 

Thus, to obtain consistency with our model (3) of quadratic 
time variation of the propagation delays T k (f), we use the 
results given in Appendix A to obtain the correspondences 
shown in Table I. 

A Delay-Rate Mismatch with Perfect Delay-Acceleration 
Match 

1 ) Derivation of Results: From (18) with a = 0, 



The corresponding r -domain impulse response h d (r, X, a) is 
obtained by inverting (22), and differs according to whether 
X^cer/2: 



i r 12 ~ 

*d(r,X,0) = — iU&(r + Xf)] dt 

1 J -T/2 

R 5 (x)dx (36) 



t r MT+(\T/2)) 
Pt*- 1 Jb.(t-{\TI2)) 



h d (r,\,a)} x>aT/2 -4 



1 



rVX* - 2oer ' 
-0, 



(30) 



otherwise 



and 



h d {T, X, a)) x< aT/2 = h d (r, X, a)) x> ar/2 + Mr, X, «) 
where 



A 2 (t,X,«H 



0, 



l -2(ir \ 8 2 / 2a 

(31) 



otherwise. 



with an equivalent expression given by (21 ) and (28) as 

**(r,A,0)=j[ j^W.(£)**(\ir)e' wtr 4r. (37) 

The frequency functions pertinent to (37) are shown in 
Fig. 3 for a block source spectrum of width B centered at f 0 , 
normalized so that 



For special cases this reduces to 
hd(r,0,0)=8(T) 



«,(r)=sinc(Br)e^<> T . 



(38) 



Arf(T-.X,0)-i 
and 



AT ,T,< ~T 



(32) 
(33) 



(.0, otherwise 

2 -or 5 



Since H d (f, X, 0) is purely real, then (refer to Fig^2) by the 
principle of stationary phase [12] the r<entroid of R d (r, X, 0) 
is the same as that of /^(for), which is at r = 0. Thus 
H d (J y X, 0) introduces no r-bias which varies with X, arid we 
need consider only the main-axis responses R a (T,0, 0) and 

£ rf (0,X,0). 
From (36) and (38), 



T^p2ar' 8 
(.0, otherwise 



<t <0 



(34) 



R d (r, 0, 0) = *,<p»= sine (fc*r) e***^, 
while from (37) and Fig. 3, 



(39) 
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TABLE 1 

Relations between Delay and Range Coefficients 



ITEM 


KINEMATIC EQUIVALENT 


MOVING SOURCE, 
FIXED RECEIVERS < s ) 


FIXED SOURCE, 
MOVING RECEIVERS 


a ,i. ' k » 1 2 


Mil ♦ r^to 








r 0k ,c 


V r 01- T 02 




<r or r o2 )/c 


A 0 " r 0l " T 02 






A o* r oi" r ie 




tr 6i " r cc >/c 
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Fig. 3. Scaled source spectrum and equivalent filter response for delay- 
rate mismatch \ and correlator integration time T. 




0 0.5 1.0 1.5 2.0 
x-f o 0 2 XT-TAMO) 

Fig. 4. Delay-rate mismatch response at r = 0 (block source spectrum, 
bandwidth B, centroid f 0 h 

*d(0,X,0) = — sinc(\Tf)df 



while at the other extreme of a low-pass block spectrum of 
video bandwidth 2? = 2/ 0 , 

Si(2nx) 



^ Si[nx(l+n)) - S/[tt*(1 - m)] 



(40) 



where 



B 



2/ 0 

N f* sin/ 
Si(x)= —dt. 

J o 1 

For a source bandwidth 5 0, the spectrum in Fig. 3 goes 
over to $(/- fo/o) and (40) becomes 

**(0, \ 0] B - 0 = sine (P 2 / 0 Xr), (41) 



*<*(0,A,0)] fl/Wo = 



27TX 



(42) 



Equations (41 ) and (42) are plotted in Fig. 4. 
We note from (40) that R d {0, X,0) is purely real. Ihus 
from Fig. 4 the 3 dB X-width (X 3 ) of either R d (0, X, 0) for a 
wide-band correlator, or of [^(OAO)! for a narrow-band 
correlator, satisfies 
= 0.88 

/0P2X3 ~ t 



for B^0 



(43) 



and is only about 1 1 percent smaller for B-+2f 0 , where / 0 is 
the centroid of the source spectrum. Thus the allowable mis- 
match (X) in delay-rate compensation for a peak 2 correlator 
loss of at most 3 dB is essentially 



IXK- 



0.44 



(44) 



To translate this tolerance into that of an allowable fre- 
quency difference, we observe that, for the component of the 
source signal at the centroid frequency /=/o, the instanta- 
neous frequency difference at t = 0 between the signals at the 
input to the multiplier in Fig. 1 is obtained from (15)-(17) as 

A/ o (0) = / o e'(0) 

~/o(&X c - A^) = / 0 ft i X (45) 

so that the abscissa of Fig. 4 may be interpreted as 

x=r-A/ o (0) (46) 

where A/ o (0) is the uncompensated frequency difference at 
/ = 0 at frequency (fe/o)- The mismatch tolerance specifica- 
tion given by (44) is then 

|A/ 0 (0)|<^ (47a) 

which, for 3 dB peak loss, corresponds to a linear phase rota- 
tion at frequency (0j/ o ) during time T of 

27rr-A/ 0 (0) = 2ir(0.44)=il58 o . (47b) 

If no delay-rate compensation is attempted, then X c = 0 in 
(17) so that 

X=-Ai/fe (48) 
and the abscissa in Fig. 4 then has a magnitude of 

x-\TfM, (49) 
which from Table I is seen to be 

*= 17/0(0, -M 
= \TfoMi(r'oi -r M )lc\ (50a) 

2 This assumes \T is small enough so that the mainlobe responses 
versus r of R d (r 9 X, 0) for the broad-band, and \R d (r 9 X, 0)\Joi the 
narrow-band correlator are unimodal with global maxima at r - 0. For 
the block source spectrum, calculations of R d (r 9 K 0) show that this 
requires that x » tnan 0 6 to 10 » 1he exact 

bound depending on B(2f 0 and the correlator type. At larger values of 
x, the mainlobe responses remain even in r but become bimodal, with 
a local minimum at t = 0. 
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for moving source and fixed receivers 

-IWw-rkJ/cl (50b) 
for fixed source and moving receivers. 

Since Pi and are within 1 percent of unity, (50a) and (50b) 
are essentially equal. 

2) Comparisons With Published Results [ 2J, fSJ, [5], { 13]: 
From (50b), Fig. 4 is seen to be equivalent to Fig. 4 of [3] , 
but with a more general abscissa. 

The approach used and the general results obtained 3 are 
identical to those given in [13] for the case where the delay- 
acceleration mismatch is zero and where real rather than com- 
plex signal notation is used. Thus the real part of (18), with 
a = 0 and &X ■+ 5, is the same as equation (3.10.2-33), p. Ill 
of [13]; and for the block source spectrum, the further nota- 
tional changes fe/o ~*fc> fa B ^ B reveal *»t our equations 
(40), (41), and (42) are identical to equations (3.10.2-35,-36, 
-38), p. 113,0/7. cit 

The equivalent filter, when no delay-rate compensation is 
used, is given by (28) and (48) and Table I 

H d (f -Ai/fe, 0) = sine (/TAi/ft) 



= sinc [(1-hlMfT] 



(51) 



which is slightly different, both from the result given in (21) 
of [5], and also from the result implied by (26b) of [2]. It 
appears that the difference from [5] stems from the definition 
of correlator clock time used there, while the difference from 
[2] rises from an approximation neglecting spectral compres- 
sion made there [p. 1 548, paragraph preceding (26a)] . 

As a final comparison we note that the "bias" in the time- 
delay estimate suggested in [2] as a consequence of ignoring 
delay rate is an artifact which vanishes if we agree that the 
correlator-delay estimate applies at the center of its data win- 
dow. The peak of R d (T y X, 0) occurs at r = 0 for any X 4 , so 
that from (17) and Table I the mean correlator-delay estimate 
is 



1 cm q * 

which for the moving source/fixed receiver scenario is 
Tcm = 'oi -r<njl c 

( ^01 ~ r 02 \ /*01 ~ ^02 \ _ 
—c—r\T-) T ° i 

s[^(-r 0 i)-r 2 (-roi)]/c 
if we neglect the acceleration terms r^rlt/lc. 



(52) 



Thus in this case the mean correlator firae-delay difference 
estimate for any X is a measure of the range-difference delay 
at a rime earlier than t = 0 by the propagation time (at t = 0) 
from the source to sensor "1." This geometric, or kinematic, 
bias is a time-lag effect arising from finite propagation veloc- 
ity, and is almost always negligible. 

Hie corresponding result for the fixed-source/moving- 



receiver scenario is 



58 (roi 



(54) 



which has a scale-factor geometric error of up to 1 percent, 
again almost always negligible. 

C Delay-Acceleration Mismatch with Perfect Delay-Rate 
Match 

1 ) Derivation of Results: From (1 8) with X = 0, we have 



i r r/2 

X d (r, 0, a) = — R s lMr + af 2 /2)] dt (55) 

1 J -Tf2 



with an equivalent expression given by (21) and (29) as 

where Z is the complex Fresnel integral (25). 

In this case both |jR d (r,0,a)| and Re [R<j(t> 0, a)] have 
their peaks at (different) nonzero values of t which vary with 



max \R d (r t 0, a)\ = l^ d (r e (a), 0, a) I 
W 

max Re[^(r, 0, a)] = Re [*„(r,(a), 0, a)] 
W 



(57) 
(58) 



so that the peak loss must be evaluated at either r e or ry, 
rather than at r = 0. These a-varying r-biases arise because the 
time shift function in the argument of the integrand in (55) is 
not odd about zero. As a consequence the equivalent filter Qn 
the sense of Fig. 2) has a nonlinear phase response over fre- 
quency, the study of which allows a good understanding of 
these effects. 
From (25) and (29), the equivalent filter is 



which we write in amplitude/phase response form as 

H d {f 0, a) = A(f 0, a) a > 
where 



(59) 



(60) 



A(f, 0, a) = j Vc*(x) + S l (*) 
< 53) . <Kf, 0, a) = tan" 1 [S(x)IC(x)] . 



(61) 
(62) 

Hie phase-delay (7$) and group-delay (T g ) functions for this 
filter are [10] 



3 We thank a reviewer foi pointing this out. 
Assuming as before that X7" is not so large that the mainlobe of 
R</(t, X, 0) becomes bimodal. 



, -»(/,0,tt) _ -«T* 



2irf 



Artx 2 



tan" 1 [S(x)/C(x)] 



(63) 
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JgU) 2n bf 



-aT 2 \<Xx) sin (tt/2) x 2 - S(x) cos (7r/2)x 2 1 
8trx [ C s (x) + 5 3 (x) " J' 



with limiting small-argument values given by 

-aT 2 



Km [?*(/)] = lim I Wl = 



24 



(64) 



(65) 



which is the average of the function (-at 2 /2) over the interval 
(-772,772). 

The amplitude, phase, and delay responses are plotted in 
Fig. 5 versus the scaled frequency variable x 2 = faT 2 /2. 

Now consider the blocksource spectrum of width B centered 
at /o. The frequency domain relations pertinent to (56) are 
similar to Fig. 3, but with sine (XT*/) replaced by H d {f, 0, a). 

For a very narrow-band signal (B^O), the compressed 
source spectrum goes over to 8(f- 0 2 / o )> (56) gives 



^,0,a)] fl . 0 =eWAr|i 

-x fish*. 



^0 



where y Q : 



(66) 



Here the envelope of the mean correlator output is 

l*d(r, 0, a)| = — VC*0> o ) + S*O> 0 ) (67) 

which is independent of r since 5 = 0, while the real part of 
the mean correlator output is 

C(yo) cos 27r/ 0 g 2 r - S(y 0 ) sin 27r/ 0 g 2 T 

yo 



Re[^(T,0,a)] = 



1 



= — Vc*(* ) ) + s>0'o) 

•COS27T/o/J a [T- ^(/o/3 2 )] 



(68) 



where ^(/ofe) is the phase delay of the equivalent filter at 
frequency (/ 0 ftt)> as given by (63) with x=y 0 . From (65), 
(66), and (68) and Fig. 5, we see that for 



J>o=/o/WV2<1.73 



(69) 



then both the envelope and the peak real response of the cor- 
relator output are attenuated by less than 3 dB, while 



aT 2 
24 



(70) 



Thus, for a very narrow-band signal with centroid frequency 
f 0 , the 3 dB tolerance value for delay-acceleration mismatch is, 
from (69), 



loci 



3.46 
MoT 2 ' 



(71) 



For the component of the source signal at its centroid fre- 
quency f Q9 the difference-frequency rate between the signals 
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Fig. 5. Amplitude, phase, phase-delay, and group-delay functions for 



- , <MM1 



at the input to the multiplier in Fig. 1 is obtained from (15)- 
(17) as 

A/ o f (0) = / o e"(0)^/o^a (72) 

so that the equivalent 3 dB tolerance for a mismatch in differ- 
ence-frequency rate is 

3.46 
T 2 



IA/ o '(0)| 



(73) 



which, at the allowable maximum for 3 dB peak loss, corre- 
sponds to a quadratic phase rotation at the edge of the cor- 
relator window of 



JrT/2 
I 
n 



(1.73) « 156° 



(74) 



where A/ 0 ' is the difference-frequency rate at frequency 0 2 / o . 

If now we allow the signal bandwidth B to increase from 
zero, then the input to the equivalent filter of Fig. 2 is given 
by (39), which has the form of a "carrier" at frequency 
C/ofo), modulated by sine (fetfr). If ji = (fl/2/ 0 ) is small 
enough that H d (f, 0, a) has essentially constant amplitude and 
linear phase over the band (fi^B) centered at (0 2 / o ), then 
[10] the output envelope is delayed by T g (p 2 fo), while the 
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output "carrier" is delayed by ^(fe/o), which is the delay 
for the peak of Re {£</(r,0, a)}. In either case, the output 
peak is approximately reduced by the "narrow-band" factor 
given by (67). 

As {i = (2?/2/ 0 ) increases further, the concepts of phase 
delay and group delay lose meaning [10] , and we must evalu- 
ate both \R d \ and Re {R d } versus t to find first the r-biases, 
and then the peak attenuations at these values of t. For the 
block spectrum this requires numerical integrations which are 
simpler if we use (38) and (55) to obtain 



,0,a)=^r f 



772 smnBp 2 (j+at 2 ll) 



-*C V?V£? ****** < 75 > 

with 

|i = Wo. 

Equation (75) was evaluated numerically for various a and a 
in the region around jc = -a/ 12 (i.e., r = -(a 7^/24)), and the 
peak delays and corresponding peak responses obtained for 
both Re{£ d (7 r ,0,a)} and U?rf(r e ,0,a)j. The results are 
shown in Fig. 6, in which the results for B = 0 are simply the 
values given in Fig. 5 for /=/ 0 ^. It is seen that the delay- 
acceleration mismatch tolerance for a 3 dB peak loss is very 
nearly independent of bandwidth for both wide-band and 
narrow-band correlators, and is essentially given by (71) or 
(73). 

We note from (72) and (73) that the abscissa in Fig. 6 may 
be expressed as 

a=A/ 0 '(0)-r 2 /2 (76) 

and that if no delay-acceleration compensation is attempted, 
then a c =0in(17) and 

a=-Aj/& (77) 



so that 



M-l/ 0 ASr 2 /2i 



(78) 



and the algebraic signs of t> and r e are reversed in Fig. 6. 
From Table I, (78) is seen to be 

\a\=\fo^(PlrZ t -ftr u m )/c\ 
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for moving source/fixed receivers (79a) 
2 

for fixed source/moving receivers (79b) 
and for 1 1 - 0* | < 0.01 , these are roughly equal. 



aM 0 /5feaT*/2 

Fig. 6. Peak shifts and losses for delay-acceleration mismatch with 
wide-band [Re(# tf )] and narrow-band l\R d \] correlators (block 
source spectrum, bandwidth B y centroid f 0 ). 

2) Comparisons with Published Results [3], [4J: For the 
narrow-band {B/2f 0 -+0) case, the peak loss for no delay- 
acceleration compensation given in (4-32) of [4] is essentially 
(68) with 

y 0 = T^ff^il = jTy/uKll (80) 

so that by (26), the sign of the phase delay r 0 (/o0 2 ) in (68) 
is reversed. In [4] the small-argument approximation (70) 
is used. 

In (44) of [3] the degradation factor for B/2f 0 -*0 is given 
as C(y 0 )ly 0 , which is erroneous since from (68) this equals 
Re[/? d (0, 0, or)] , which is neither the real output peak nor the 
peak of the envelope of the output. This error appears to have 
arisen because the term involving the Fresnel sine integral was 
ignored in passing from (43) to (44), op. cit. 

These results are plotted in Fig. 7 for comparisons. 

IV. Summary 

We have considered the problem of how, and how well, to 
compensate for the effects of source or receiver motions on 
the output of a cross correlator used to estimate delay differ- 
ence across a baseline. For a quadratic delay difference A(r) 
during the correlator integration period T given by 

A(r) = A 0 + tit + Ai f r 2 /2, 1 1\ < Til, (4) 

the required correlator compensation was found to be essen- 
tially a quadratic delay modulation A c (r) inserted in one input 
channel before multiplication 



A c (0 = r c + \ c f + a c r 3 /2, |r| < 7/2. 



(6) 



The approach used was to assume a stochastic source signal 
and to determine the mean (ensemble average) correlator out- 
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Fig. 7. Comparison of results for peak output loss of correlator with 
quadratic-delay mismatch and negligible signal bandwidth. 



(*-•)• 



put as a function of the mismatches 
X=X C - Ai/ft 

ot = a c - A2/02, 0 2 =1-Tot 



(17) 



in compensating for delay -difference rate and acceleration. 
The use of analytic signal representation allowed broad-band 
and narrow-band correlators to be treated simultaneously. 

For no more than 3 dB reduction in the peak of the mean 
correlator output with a narrow-band signal and a perfect 
delay -acceleration match (a= 0), the delay-rate mismatch (X) 
must satisfy 



I X21 < or irA/ 0 (0)|<0.44, 



(44,47) 



while with a perfect delay-rate match (X = 0), the delay- 
acceleration mismatch (a) must satisfy 



I a 7^/2 1< 



1.73 
Pifo 



or 



«T»2 

— A/ 0 '(p) 



< 1.73 



(71,73) 



where / 0 is the centroid frequency of the source spectrum, T 
is the correlator integration time, and A/" o (0), A/ 0 '(0) are the 
frequency difference and frequency-difference rate at fre- 
quency ft/ 0 and f=0, at the input to the correlator multi- 
plier. For broad-band signals, these mismatch tolerances be- 
come about 1 1 percent tighter. 

For given values of X and a, these mismatch tolerances pre- 
scribe the maximum allowable correlator integration time 
r max as the smaller of 



0.44 



and 



I 3.46 



(t-T) 




<t) SOURCE 



~rJ*0) RECEIVER 

(t- T 



Fig. 8. Assumed propagation geometry. 

Appendix A 
Signal Range-Delay Variation 
with Moving Source or Receiver 
The aim here is to describe the time-varying signal delay 
r = r(0 in terms of the time-varying range r(r) when the 
source or receiver, or both, move. Assuming direct-path prop- 
agation in a horizontal plane in an idealized medium with 
constant speed of propagation c, the geometry is indicated in 
Fig. 8. 

The source transmits a waveform s(t- r) at time (t- t) 
when the range is r(t - t). Neglecting attenuation and distor- 
tion in the medium, this point value of the source waveform 
travels along the path cr(r) and impinges on the receiver at 
time t when the range is r(f), so that the received waveform at 
time t is 

s ,(r) = s(r-r(0). (AD 

For a fixed source and moving receiver (FS, MR) then d s = 0 
in Fig. 8, and 

cr(t)-r(0. (A2) 
For a moving source and fixed receiver (MS, FR) then d r = 0 
so that 

cr(0 = r(f-T(f)). (A3) 

When both source and receiver move (MS, MR), then by the 
law of cosines 

cr(0 = r(f) Vl - 2d s cos 0,/r(r) + te/K')) 2 (A4) 

where d s =d s (t,T) is the straight line connecting the source 
positions at times (r - t) and (f), and 0 S = d 5 (t, r) is the angle 
as indicated between d s and r(f). With an instantaneous 
source velocity schedule given by »„(f), v sy (t) during (t - r, r)> 
then the average source speed over this interval is 

V s =>/V*+V* = V 8 (! 9 t) (A5) 

with 



1 v 



:«)««= *i(f,T) 



If the correlator does not attempt to compensate for delay- 
difference accelerations, then ct c = 0 and a is determined by 
the kinematics Ao and may be large, leading to a small value 
of ^mox- A similar conclusion holds for a correlator which 
does not compensate for delay-difference rates. 



and similarly for V y . The length d s is then 

d a =V s T(t) 
where 

r(0 = — 



(A6) 



(A7) 



(A8) 
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so that 



where from (A11),(A15), ■ • • ,(A17), 



d V 

0.01 



(A9) 



To — 



r(0 c 

and a first-order expansion of (A4) yields 

cr(f) a r(r) - (V s cos $ 3 ) r(r). (Al 0) 

Equations (A2), (A3), and (A10) relate r(r) to r(t) for the 
(FS, MR), (MS, FR), and (MS, MR) cases. We now concentrate 
on the first two cases on the assumption that the range r(t) 
varies quadratically with time over the observation interval 



1 +rj>/c 



To = 



rile 



r'olc 



(A19) 



r(t) = r 0 +r' 0 t + ^t\ |<|<rJ2. 
1) FS, MR-From (A2) and (Al 1), 

r(f)=^ = ro+ri/ + r;-y 
c * 



where 



r 0 

To = 7 



r > -IS. 



(All) 



(A12) 



t5 



2) MS, FR— From (A3) and (All), 



(A13) 



(A14) 



and where we have assumed that 
ro-^o=A-To) = '\0)=>o. 

Equations (A13) and (A19) relate the coefficients of the 
quadratic signal-delay function to the coefficients of the 
quadratic range function for the two scenarios (FS, MR) and 
(MS, FR). Note that if r" = 0 all approximations become 
exact for the (MS, FR) case. 

Appendix B 

Delay-Difference Compensation Approximations 
With r k (t) and A(f) - T,(f) - r 2 (0 assumed quadratic in time 
(3), (4), then for an arbitrary delay-compensation modulation 
A c (t), the delay-difference function at the input to the multi- 
plier in Fig. 1 is the negative of e(r), as given by (15). Com- 
bining (3), (4), and (1 5) gives 

6(0= ~ [1 " ria ~ r«r] A e (f)- A(0- (Bl) 

For perfect compensation we require e(r) = 0 over (-7/2, 
F/2), so that A c (r) is the root of a quadratic 



where we temporarily suppress the explicit time-function 
notation. The solution of (A14) for r yields 

since the second term under the radical in (A 15) is of the 
order of c' 1 times the change in r' during t, and is thus much 
less than 1. 
Successive differentiation of (A14) yields 



A c (0 = 



t m 



(B2) 



In nearly all practical cases the second term under the radical 
ta(B2)is«l,sothat 



Ac(0 



a(o „ m . t 

0a - r'nt fa 



(B3) 



t n 

r - rr 



■4 (>-•) 



r 



(A16) 



(A17) 



The signal delay function for the (MS, FR) case is then ap- 
proximately quadratic for a quadratic range variation 



T(r)-ro+r 0 r + y r 5 



(A18) 



therefore, if A(f) is quadratic in time, then the optimum A e (f) 
is very nearly quadratic in time. 5 Thus we assume that A e (r) 
has the quadratic form given by (6), and when (3), (4), (6), 
(IS) are combined, we obtain 

e{t) = A +Bt+2 f 4 +Et i +Fr 4 (B4) 

where 

B = h\ c ~K- [tm(1- A c )r c ] 

D = faa c - K ~ [rm {X c (2 - X c ) - r c a c }) 

£=[-TS,(l-X c )a c /2] 

F«[r&«3/8]. (B5) 
5 This is exact if t'qz = 0. 
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By neglecting the terms in brackets in (B5), we obtain (16). 
In any particular scenario the effect of these approximations 
can be estimated by expanding R s [€(t)] to low order in e(t) 
and integrating (14). In nearly all cases the approximations 
leading from (B5) to (16) are very good. 
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